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Abstract 

In connection with recent and proposed experiments, and new theoretical results, my 
previous calculations of the Lamb shift in muonic hydrogen will be reviewed and compared 
with other work. In addition, numerical results for muonic deuterium and helium will be 
presented. Some previously neglected (but very small) effects are included. 

Introduction 

The energy levels of muonic atoms are very sensitive to effects of quantum electrodynamics 
(QED), nuclear structure, and recoil, since the muon is about 206 times heavier than the 
electron [1]. 

A recent measurement of the Lamb shift in muonic hydrogen [2] has stimulated great 
renewed interest in the energy levels of muonic atoms. A number of theoretical analyses 
of the Lamb shift (the 2p-2s transition) in light muonic atoms have been published [31 IH 
El El [3 m m Uni im HSJ before the experiment was performed. The present paper repeats 
the independent recalculation of some of the most important effects [6] for hydrogen and 
deuterium [13], including some previously neglected (but very small) effects and extends 
the numerical calculations to the case of muonic helium. Some recent results by other 
authors are also discussed. 

In the numerical calculations the fundamental constants from the most recent 
CODATA compilation ([E]) are used (with the exception of the fine structure constant, 
which was revised more recently), i.e.: 

he, m^, me, m„ = 137.0359991, 197.32696 MeV-fm, 105.658367 MeV, 0.5109989 MeV, 
931.494030 MeV, respectively 

The changes in these constants compared with CODATA 2002 ([I5]) are too small to 
make any relevant difference in the results, but the most recent values are used anyway. 

Also, the following properties of the proton and deuteron were used: 
mp = 938.27203 MeV/c^, /Xp = 2.79285 /x^, and i?p = 0.875(7) fm (from [15]); this agrees 
with a recent determination by the group at Jefferson Laboratory [TH]. Other recent values 
have been reported, for example 0.8418(7) fm [2] and 0.877(8) fm [T7]. In some cases calcu- 
lations were performed for several values of Rp. In addition (see [H]) = 1875.613 MeV/ c^, 
lid = 0.85744 /iTv, i?d = 2.139(3) fm (from [I5]), or i?^ = 2.130(3) fm, ([E]). For the he- 
lium isotopes m„ = 3727.379 Me V/c2, i?„ = 1-676(8) fm ([19]), or i?^ = 1.681(4) fm ([20]); 
"^He3 = 2808.391 MeV/c^ and /ijifes = — 2.1275//^^. An old measurement of i^j^es = 1-844 fm 
has been superseded by newer results which range from 1.959(34) fm [22] or 1.961 fm 
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to 1.975(4) fm (other results are cited in these papers). Here an average value of 
1.966(10) fm will be used for an estimate, since there is still some uncertainty in the 
calculations used in the newer determinations. In some cases, older values for the proton 
and deuteron radii were used to make comparison with earlier results (O [13] easier. 

The deuteron has spin 1 and thus has both magnetic and quadrupole moments. The 
quadrupole moment of the deuteron is taken to be Q = 0.2860(15) fm^ [23 EH EZ]- A 
newer result |28] is nearly the same, but more precise. 



Vacuum Polarization 

The most important QED effect for muonic atoms is the virtual production and anni- 
hilation of a single e'^e~ pair. It has as a consequence an effective interaction of order 
aZa which is usually called the Uehling potential ([221 ED]- This interaction describes 
the most important modification of Coulomb's law. Numerically it is so important that it 
should not be treated using perturbation theory; instead the Uehling potential should be 
added to the nuclear electrostatic potential before solving the Dirac equation. However, 
a perturbative treatment is also useful in the case of very light atoms, such as hydrogen. 

Unlike some other authors, we prefer to use relativistic (Dirac) wave functions to 
describe the muonic orbit. This is more exact, and as will be seen below, it makes a 
difference at least for the most important contributions. Relativistic recoil corrections 
to this treatment will be calculated according to the prescription given in [T]. The wave 
functions are given in the book of Akhiezer and Berestetskii [31j and will not be given 
here. In perturbation theory, the energy shift due to an effective potential AV is given 
by 

-| poo poo 

A^n. = ^ - / q'dqAViq)- / drUqr)[Fl + GU (1) 
Jo Jo 

where and Gnn are the small and large components of the wave function, n is the 
principle quantum number and k is equal to — £ — 1 if j = £ + | and +£ if j = C — \. 
AV{q) is the Fourier transform of the physical potential. 

/•oo 

AV{q) = An- • jo(gr) ■ AV{r) dr (2) 

Jo 

1 f°° 

^V{r) = q'-Uqr)-AV{q)dq (3) 

Jo 

As is well-known [T], the Uehling potential in momentum space is given by 

Vuehiiq) = -l^Ml . GM ■ F{<P) = -An{aZ) • G^(g) ■ UM 
where Ge is the proton charge form factor, sinh(0) = g/(2me) and 

F{<P) = ^ + (coth2(0) - 3) • [1 - ■ coth(0)] (4) 

U2{q) is also defined in [T]. If the correction to the transition 2pi/2 — 2si/2 is calculated 
in lowest order perturbation theory using nonrelativistic point Coulomb wave functions, 
the result is 205.0074 meV, in agreement with other authors^. 
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The same procedure was used to calculate the two-loop corrections; the corresponding 
diagrams were first calculated by Kallen and Sabry [32]. The Fourier transform of the 
corresponding potential is given in [HIT]. The result for a point nucleus is 1.5080 meV. 

In momentum space including the effect of nuclear size on the Uehling potential is 
trivial, since the corresponding expression for AV{q) is simply multiplied by the form 
factor. The numbers obtained were the same for a dipole form factor and for a Gaussian 
form factor, provided the parameters were adjusted to reproduce the assumed rms radius 
of the proton. The correction can be regarded as taking into account the effect of finite 
nuclear size on the virtual electron-positron pair in the loop. The contribution of the 
Uehling potential to the 2p-2s transition is reduced by 0.0079 meV with a proton radius 
of 0.842 fm [2], and by 0.0082 meV with a proton radius of 0.875 fm [15]. This result is 
consistent with the number given in [3] (eq.(266)). The contribution due to the Kallen- 
Sabry potential is reduced by 0.00007 meV. 

These numbers are modified somewhat when relativistic (Dirac) muon wave functions 
are used. The numerical values given below were calculated as the expectation value of 
the Uehling potential using point-Coulomb Dirac wave functions with reduced mass: 





point nucleus 


i?p=0.875fm 


Uehling 
Kaellen-Sabry 


2pi/2 - 2Si/2 2^3/2 - 2Si/2 

205.0282 205.0332 
1.50814 1.50818 


2pi/2 - 2Si/2 2p3/2 - 2Si/2 

205.0198 205.0248 
1.50807 1.50811 



The effect of finite proton size calculated here was checked for several values of the 
proton radius, ranging from 0.842 fm to 0.892 fm. It can be parametrized as -O.OllO(r^). 
The contribution due to two and three iterations has been calculated by |H |35] and 
|34] . respectively, giving a total of 0.1507meV. An additional higher iteration including 
finite size and vacuum polarization is given in ref. |1] (equations(66) and (67)) and ref. [3] 
(equations(264) and (268)). This amounts to -0.0164(r^), and is verified in Appendix B. 
This would mean that the finite-size contributions to vacuum polarization in muonic 
hydrogen can be parametrized as — O.OllO(r^) — 0.0164(r^), giving a total of — 0.0274(r^). 
Numerically this is -0.0209(6) meV if the proton radius is 0.875 fm, or -0.0194 meV if the 
proton radius is 0.842 fm. 

Up to now the higher order effect including nuclear size has only been calculated for 
hydrogen. The contribution for other light nuclei is given in Appendix B. 

Corresponding numbers for muonic deuterium, calculated as the expectation value of 
the Uehling potential using point-Coulomb Dirac wave functions with reduced mass are: 





point nucleus 


i?d=2.139fm 


Uehling 
Kaellen-Sabry 


2pi/2 - 2Si/2 2p3/2 - 2Si/2 

227.6577 227.6635 
1.66622 1.66626 


2pi/2 - 2Si/2 2p3/2 - 2Si/2 

227.5986 227.6043 
1.66577 1.66581 



The effect of finite nuclear size calculated here can be parametrized as — 0.0130(r^). 
However higher iterations will change these results. The correction for a point nucleus has 
been calculated for muonic deuterium [22], giving an additional correction of 0.1718 meV. 
The effect of finite size as described in refs. [21 S] is given in Appendix B and results 
in an additional contribution to the binding energy of the 2s-state of — 0.02062 (r^), or 
— 0.0936 meV for a deuteron radius of 2.130fm. 
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Corresponding numbers for muonic ^He and ^He, calculated as the expectation value 
of the Uehhng potential using point-Coulomb Dirac wave functions with reduced mass 
are: 





point nucleus 




^He 

Uehling 
Kaellen-Sabry 


2pi/2 - 2Si/2 2^3/2 - 2Si/2 

1666.305 1666.580 
11.5731 11.5752 


2pi/2 - 2Si/2 2^3/2 - 2Si/2 

1665.381 1665.656 
11.5658 11.5680 


3He 

Uehling 
Kaellen-Sabry 


2pi/2 - 2Si/2 2^3/2 - 2Si/2 

1642.412 1642.682 
11.4107 11.4128 


2pi/2 - 2Si/2 2^3/2 - 2Si/2 

1641.337 1641.607 
11.4024 11.4045 



Here the nuclear radii were taken to be -R//e3=l-966 fm and i?//e4=l-681 fm. The ef- 
fect of finite nuclear size calculated here (and with other radii) can be parametrized 
as -0.3297(r2) for ^He and -0.3176 (r^) (-0.3151 (r^) if the two-loop contribution is 
neglected) for ^He. Higher iterations change these results. In pQ they were calculated 
to be about 1.70 meV, including the effect of finite nuclear size, for '^He and 1.4 meV ^He. 
A more recent calculation of the correction for a point nucleus has been made for muonic 
^He [35], giving a correction of 1.709 meV. The effect of finite nuclear size on the shift of 
the 2s-state is given in Appendix B (for both isotopes). 

The mixed muon-electron vacuum polarization correction was recalculated and gave 
the same result as obtained previously, namely 0.00007 meV. [MIE]- For deuterium, the 
contribution is 0.00008 meV, for ^He it is 0.00200 meV, and for ^He it is 0.00208 meV. 
For the helium isotopes, neglecting finite nuclear size would increase the contribution by 
0.0001 meV. 

The Wichmann-KroU [23] contribution was calculated using the parametrization for 
the potential given in [T]. The result obtained for hydrogen is -0.00103 meV, consistent 
with that given in [3]. For deuterium, the contribution is -0.00111 meV. Values for both 
helium isotopes have also been calculated. The contribution is -0.01984 meV for ^He and 
-0.01969 meV for ^He. The reason for the negative sign is discussed in the review of Eides 
et al. 0. 

The sixth order vacuum polarization corrections to the Lamb shift in muonic hydrogen 
have been calculated by Kinoshita and Nio [51]. Their result for the 2p-2s transition (in 
hydrogen) is 

AE^^) ^ 0.00761 meV 

Recently Karshenboim et al. [351 EZ] have recalculated all of the relevant graphs, and 
reproduced these results for hydrogen, including a correction that reduced the contribution 
to 0.00752 meV. In addition they performed a recalculation of the virtual Delbriick effect 
that is consistent with, but probably more accurate, than the results given in [38] and 
[1] , possibly due to a more accurate numerical evaluation of the sevenfold integrals over 
Feynman parameters. The published results are given as the sum of this, the Wichmann- 
KroU contribution and the previously uncalculated light by light contribution. In in the 
case of muonic hydrogen and helium, the virtual Delbriick contribution very nearly cancels 
the Wichmann-Kroll contribution, so the total result can be expected to be small. 
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Numerically, they found the following total " light-by- light" corrections (i.e. the sum 
of the Wichmann-Kroll, virtual Delbriick and previously uncalculated term), in meV: 



AE{fiH) -0.00089(2) 
AeIiiD) -0.00096(2) 
AE{fi^He) -0.00136(6) 

The hadronic vacuum polarization contribution has been estimated by a number of 
authors |39lll0l|3]. It amounts to about 0.012 meV in hydrogen, 0.013 meV in deuterium, 
0.219 meV in ^He and 0.225 meV in ^He. An uncertainty of about 5% should be expected. 



Finite nuclear size and nuclear polarization 

The main contribution due to finite nuclear size has been given analytically to order (aZ)^ 
by Friar [H]. The main result is 



2aZ ( aZnir 



n 



aZrrir 



{r^)^2) + {aZf{FREL + mlF, 



NR 



(5) 



Radiative corrections to the main term have been calculated by Fides and Grotch [12] . 
They contribute an additional correction of 



AEr, 



2aZ ( aZrUr 



n 



a^Z{2'i/A - 41n(2) - 3/4) ■ (r^) 



This correction (of order a{aZ)^) amounts to an additional (fractional) correction to the 
main term equal to 2.2275a^Z(r^) = 1.1862 x 10~^Z(r^). For muonic hydrogen, the main 
coefficient of (r^) is -5.1973 meV fm^^ without the radiative correction; it is modified to 
-5.1979 meV fm~^ with the correction. For the main term, the shift is -3.979±0.076meV if 
the proton rms radius is 0.875±0.007) fm and -3.685±0.008 meV if the proton rms radius 
is 0.842±0.001) fm. The radiative correction increases these numbers by 0.0005 meV. 
Assuming a proton rms radius of 0.875 fm, the correction from the first two terms in 
Eq.([5]) to the 2s-level is 3.956 meV for a dipole form factor and 3.958 meV for a Gaussian 
form factor. The parameters were fitted to the proton rms radius. The difference is due 
to the second term in Eq.dS]). It contributes 0.009126 (r^)(2)meV (or 0.0232 meV for a 
dipole form factor and 0.0212 meV for a Gaussian form factor). 

Pachucki [5] has estimated a correction similar to the second term (proportional to 
(r'^)(2)) in Fq.©. Friar and Sick [13] have shown that the results are equivalent, at least 
for static potentials. 
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Since the second term in Eq.([S]) is numerically important, it has been investigated 
further, at least for the case of hydrogen. The third Zemach moment was calculated in a 
model-independent manner from electron-proton scattering data by Friar and Sick 



with the result {r^)(2) = 2.71(13)fm^, for a contribution of 0.0247 meV, with an uncertainty 
of 0.0012 meV. and more recently, using new experimental results, by Distler et al. |44] . 
with the result {r^){2) = 2.85(8)fm^, for a contribution of 0.0260 meV, with an uncertainty 
of 0.0007 me V. The relationship between the third Zemach moment and the charge radius 
will be discussed below. 

Since the coefficient of (r^) is important for the determination of the nuclear radius, 
the terms of order (aZ)^ in Eq.(|5]) that depend only on this quantity are separated from 
the other terms. This affects only [aZYF^EL- Details will be given in Appendix B. The 
result is an extra contribution to the energy shift of the 2s-state equal to 

2 {aZfml, 2, . 35 , , 

(7 is Euler's constant). This results in an energy shift of — 0.00181 (r^) = — 0.00138 meV 
in hydrogen. The value of the coefficient is compatible with the estimated contribution 
of — 0.0016(r^) given in ref. [2j. The other remaining terms (of order {aZY) given in 
[H] contribute 0.000123 meV if an exponential charge distribution is used. Some details 
are given in Appendix B. This estimate includes all of the terms in Eq.([5]), while other 
authors [3l [5] give only some of them. 

Friar [H] also found finite size contributions to the binding energy of the 2p-levels, , 
amounting to 



{aZY ( aZm 



AKpi/2 = ■ (1 - l/n') 



n 



~2~ ^ 15 



(6a) 



, , faZmrY (^2 - l){r^) 
Ai?„,3/2 = {aZ) (^— -^J ■ ^ (6b) 

If terms of order {aZY are included in the coefficient of (r^) are to be taken into 
account consistently, the contribution to the 2pi/2 energy level must also be included. 
For hydrogen, this is -0.0000519(r^) = -4-10^^ meV. Subtracting this contribution from 
the order {aZ)^ shift of the 2s level results in a contribution to the 2si/2-2pi/2 transition 
of 0.00131 meV, in good agreement with previous estimates There is no contribution 
to the coefficient of (r^) for transitions to the 2p3/2 state. The fine structure is, of course, 
affected. The contribution to the 2pi/2 energy level of deuterium, gives a correction of 
-0.0000606(r2) =-2.8-10-^ meV. 

The contribution to the energy of the 2p3/2 state is much smaller, and is proportional 
to (r^); the relation to the rms-radius of the nucleus is model dependent. Analytic 
expressions for (r^) in terms of (r^) can be derived from the table given in ref. [H]. For an 
exponential charge distribution, (r^) = 5((r^))^/2 and for a Gaussian charge distribution 
(r^) = 5((r^))^/3. For the helium isotopes these corrections are significant for both levels. 

These contributions, with numerical values for the coefficients and energy shifts for 
deuterium and helium will be given in Appendix B. 

For the extraction of a single radius parameter from experimental data, it is useful, 
and has become customary, to write the Lamb shift in the form 
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It is straightforward to determine B. There are several contributions, and these will be 
given in Appendix B, with numerical values for all cases. However, the value of C depends 
on a model for the charge distribution, since it depends on the ratio 

For Gaussian and exponential charge distributions, / can be calculated analytically |41j . 
with the result / = 32/(3 v^Svr) = 3.4745 for a Gaussian charge distribution, and 
/ = 105/(16a/3) = 3.7889 for an exponential charge distribution, which corresponds to 
a dipole form factor. The corresponding coefficients of ((r^))'^/^ (in units meVfm"^) 
are, for hydrogen, 0.03458 (dipole form factor) and 0.03171 (Gaussian form factor). 
When the third Zemach radius is determined from experiment, the value of / must 
also be determined from the same experimental data. Using the results of Friar and 
Sick HHIlS] gives / = 3.780(52), while the newest evaluation P] gives / = 4.18(13), The 
corresponding coefficients of ((r^))^/^ are 0.0345(5) and 0.0383(12), respectively. In the 
newest CODATA compilation [T3] a value of / = 3.4(2) is given. This is almost certainly 
too small in view of the previous values discussed. An average value of the two model 
independent determinations corresponds to / = 4.0(2) and would result in a coefficient 
of 0.0365(18) meVfm"^. The uncertainty in the precise value of the coefficient will have 
consequences for the theoretical uncertainty in the value of the proton radius, evaluated 
as in the recent experiment [2]. The value of A also has a theoretical uncertainty, which 
may have been slightly underestimated, (see the summary tables and Appendix C). 

The recommended contribution due to the third Zemach moment (expressed in terms 
of ((r2))3/2) for deuterium is (0.0112/ = 0.0448(22)): 0.0448(22) meVfm-3((r2))3/2. 

For the helium isotopes a Gaussian charge distribution is a fairly good approximation 
IH], so a value of use / = 3.5(.l) should be a sufficiently good approximation. Then 
C = 1.40(4) for ^He and 1.35(4) for ^He. 

The contribution due to nuclear polarization (in hydrogen) has been calculated by 
Rosenfelder [46] to be 0.017 ± 0.004 meV, and by Pachucki [5] to be 0.012 ± 0.002 meV. 
Other calculations [171 SB] give intermediate values (0.013 meV and 0.016 meV, respec- 
tively). A very recent calculation [IH] gives a unified treatment of elastic and inelastic 
contributions. The total given there is very close to the total that would be obtained in this 
work if only the inelastic contribution (0.0127(5)) is included in the nuclear polarization 
correction, and the Zemach moment contribution is interpreted as the sum of the other two 
contributions given there. However, the estimated error in the polarizability correction 
has been enlarged to be comparable to the total uncertainty. 

The contribution due to nuclear polarization in deuterium has been calculated by 
Leidemann and Rosenfelder [50] to be 1.50 ± 0.025 meV (see also [51]). A newer calculation 
[52] gives 1.680 ± 0.016 meV. The author also claims that some inelastic contributions to 
the polarizability cancel the contribution due to the Zemach moment. 

An estimate of the nuclear polarization contribution for the helium isotopes is given 
in [H IS3]- It amounts to 3.1 ± 0.6meV for ^He and 4.9 ± 1.0 meV for ^He. 

Relativistic Recoil 

As is well-known, the center-of-mass motion can be separated exactly from the relative 
motion only in the nonrelativistic limit. Relativistic corrections have been studied by 
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many authors, and will not be reviewed here. The relativistic recoil corrections summa- 
rized in [1] include the effect of finite nuclear size to leading order in 171^/171^ properly. In 
fact the effect calculated here is probably accurate to order aZm^/niN 

Up to 2004 this method was not used to treat recoil corrections to vacuum polarization. 
The recoil correction to the energy shift due to the Uehling potential can be included 
explicitly, as a perturbation correction to point-Coulomb values, in a manner very similar 
to that described in P (and in [2]). This was first described in [6]; the result in the 
case of hydrogen for the total relativistic correction agrees very well with an independent 
calculation based on a generalized Breit equation [56]. The basic formulas and numerical 
results are given here. Details are given in Appendix A. To leading order in l/m]sr, the 
energy levels including a shift beyond that due to using the reduced mass in the Dirac 
equation are given by 



(7) 



where Er is the energy level calculated using the Dirac equation with reduced mass and 
Bo is the unperturbed binding energy. Also 

h{r) = -P,{r){P,{r) + -Q,{r)) - lQ,{r)[P,{r) + ^Q,{r)] (8) 



Here 



/■oo 

Pi(r) = AixaZ / r' p{r')dr' 

J r 



Q2{r) = AnaZ / r p(r )dr 
Jo 

Q^(r) = AnaZ / r'^p{r')dr' 
Jo 



■ - V{r) - rV'{r) 
:rV(r) 



:rV(r) - 2rV(r) + 6 / r"V{r')dr' 

'0 



(9) 



Details for the case of vacuum polarization are given in Appendix A and in Ref. 

To obtain the full relativistic and recoil corrections, one must add the difference 
between the expectation values of the Uehling potential calculated with relativistic and 
nonrelativistic wave functions, giving a total correction of 0.0169 meV for muonic hydro- 
gen. This is in very good agreement with the correction of .0168 meV calculated by 





Main term 


Correction 


Total 


hydrogen 


0.0210 


-0.0041 


0.0169 


deuterium 


0.0229 


-0.00246 


0.0204 


3He 


0.527 


-0.032 


0.495 


^He 


0.534 


-0.0248 


0.509 



Table 1: Total relativistic corrections due to vacuum polarization to the Lamb 
shift in light muonic atoms (in meV). The main term is the difference 
between the vacuum polarization contribution calculated with relativistic and 
nonrelativistic wave functions. The correction is calculated here. 
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Veitia and Pachucki [55]. The treatment presented here has the advantage of treating the 
main contribution relativistically and applying a small correction that can be calculated 
using first order perturbation theory. The Kallen-Sabry contribution to the main term 
was included. For hydrogen, this adds 0.0002 meV to the main term. The results for all 
cases are summarized in Tabled] 

A recent paper by Jentschura \E7\ overlooks the fact that a correction due to relativistic 
recoil was included in previous work [6J and adds an additional "correction" of -0.0039 meV 
to the Lamb shift of muonic hydrogen in Eq. (3.3). This amounts to double counting this 
contribution, which is listed in Table|3] (and in ref [6]) as a "recoil correction to VP". 

This method is also used to calculate relativistic recoil corrections to the shift due to 
finite nuclear size when a model for the charge distribution is given. This was done by 
Friar [41J; The recoil correction is given by: 

= ( j r)(2)d)^o- 

Numerically the contribution due to finite nuclear size to the recoil correction for the bind- 
ing energy of the 2s-level is -0.013 meV in hydrogen, and -0.019 meV in deuterium. The 
factor l/rriN is replaced by l/(m^ + 771^), also consistent with the calculations presented 
in IH]. In the case of hydrogen, the contribution is approximately 0.004 meV smaller if a 
Gaussian charge distribution is used instead of an exponential charge distribution. With 
the exponential charge distribution, the value sf the contribution ranges from 0.0131 meV 
to 0.0139 meV as the charge radius varies from 0.842 fm to 0.895 fm. A calculation with 
the value given by Distler et al. [H] is also in this range. This contribution is thus nearly 
constant for hydrogen. For muonic ^He the recoil contribution is -0.267 meV for a Gaussian 
charge distribution and an rms-radius of 1.681 fm. Varying the radius by ± 0.02 fm changes 
the contribution by ±0.003 meV. For muonic ^He this recoil contribution is -0.404 meV 
for a Gaussian charge distribution and an rms-radius of 1.966 fm. 

The review by Eides et.al [3] gives a better version of the two photon recoil (Eq. 136) 
than was available for the review by Borie and G. Rinker [1]. Evaluating this expression for 
muonic hydrogen gives a contribution of -0.04497 meV to the 2p-2s transition in hydrogen, 
-0.02656 meV in deuterium, -0.4330 meV in ^He, and -0.5581 meV in ^He. 

Higher order radiative recoil corrections give additional contributions These 
correspond to (their) Table 9 and two additional terms in Table 8. Numerically, the 
contributions are -0.01003 for hydrogen, -0.00302 for deuterium, -0.04737 for *^He and 
-0.08102 for ^He. 

An additional recoil correction for states with i ^ has been given by [58] (see also 
[3]). It is 

Note that l/(j + 1/2) - l/(£ + 1/2) = l/(/t(2£ + 1)). When evaluated for the 2p-states of 
muonic hydrogen, one finds a contribution to the 2p-2s transition energy of 0.05748 meV 
for the pi/2 state and -0.02873 meV for the P3/2 state. For ^He the contribution is 
0.12654 meV for the pi/2 state and -0.06327 meV for the P3/2 state. For integer spin 
nuclei, there is an additional shift of the 2s-state, amounting to —aZ{aZmr/nY/{2mj^) 
[55] . Numerical values for the energy shift of the 2pi/2 and 2p3/2 states of deuterium 
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are: 0.01681 meV and -0.00840 meV, respectively. The total contribution for the 2pi/2-2s 
transition is then 0.06724 meV; the fine structure is unaffected. Numerical values for ^He 
are 0.07379 meV for the 2pi/2 state, -0.03690 meV for the 2p3/2 state, and -0.22137 meV 
for the 2s state. 

A final point about recoil corrections is that in the case of light muonic atoms, the mass 
ratio m^/mjsi is considerably larger than the usual perturbation expansion parameter aZ. 



Muon Lamb Shift 

For the calculation of muon self-energy and vacuum polarization, the lowest order (one- 
loop approximation) contribution is well-known, at least in perturbation theory. For a 
point nucleus, and neglecting possible contributions due to vacuum polarization, so that 
V'^V = —4:TTaZp (p is the nuclear charge density) is approximated by a delta function, 
giving V = —aZ/r and 

(VV) = -AnaZ\i;^sm'Sm = -4Za(^^)'5,o 

n 

one finds a contribution of -0.66788 meV for the 2si/2 — 2pi/2 transition and -0.65031 meV 
for the 2si/2 — 2^3/2 transition. For deuterium, the corresponding contributions are given 
by -0.77462 meV for the 2si/2 — 2pi/2 transition and -0.75512 meV for the 2si/2 — 2^3/2 
transition. These numbers include muon vacuum polarization (0.016844 meV for hydro- 
gen, 0.019682 meV in deuterium, 0.33225 in ^He, and 0.34132 in ^He), and an extra term 
of order [Za)^ as given in [3] for the 2s-state: 

which contributes -0.00443 meV for hydrogen and -0.00518 meV for deuterium. For the 
helium isotopes the contributions are -0.17490 meV in ^He and -0.17967meV in ^He. These 
results, and the higher order corrections [H [36] are summarized in TableEJ 
The higher order contributions can be written in the form 

where -F2(0) = a^; the higher order contributions (fourth and sixth) can be taken from 
the well-known theory of the muon's anomalous magnetic moment: 

^2(0) = = a/27r + 0.7658(a/7r)2 + 24.05(a/7r)3. The fourth order contribution to F{(0) 
is 0.46994(a/7r)2 + 2.21656(a/7r)2 = 2.68650(a/7r)2 [H ES]. These expressions include 
the fourth order electron loops [_54J which dominate the fourth order contribution. The 
contribution of the electron loops alone is -0.00168 meV for the 2si/2 — 2pi/2 transition 
and -0.00159 meV for the 2si/2 — 2^3/2 transition. Including the rest, which is the same as 
for the electron [T], gives -0.00169 meV and -0.00164 meV, respectively. Since it is no more 
difficult to include the complete contribution to F{{0) and -^2(0) in the analysis, there is 
no reason not to do so. An additional contribution due to including the Kallen-Sabry 
correction with muon loops effectively adds — (41/162)(a/7r)^ to F[{0) [3jj. This has been 
included in the higher order corrections. 
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Transition 


2pi/2 - 2Si/2 


2^3/2 - 2Si/2 


Hydrogen 
second order 
higher orders 


-0.667882 
-0.001714 


-0.650313 
-0.001647 


Total 


-0.669596 


-0.651960 


Deuterium 
second order 
higher orders 


-0.774616 
-0.002001 


-0.755125 
-0.001926 


Total 


-0.776617 


-0.757051 


3He 

second order 
higher orders 


-10.827368 
-0.033749 


-10.504167 
-0.032515 


Total 


-10.861117 


-10.536672 


^He 

second order 
higher orders 


-11.105708 
-0.034663 


-10.776650 
-0.033406 


Total 


-11.140370 


-10.810056 



Table 2: Contributions to the muon Lamb shift {E{2pi/2) — ^(2si/2)) iu muonic 
hydrogen, deuterium, ^He, and ^He, in meV. 



The numbers given in Table[2]were calculated assuming that V'^V can be approximated 
by a delta function. However, the potential should be corrected for vacuum polarization 
due to electron loops, and, at least for s-states, for the effect of finite nuclear size. This has 
been done a long time ago (0 E51 ESI EZ]) for the hyperfine structure, but up to now not 
completely for the " muon Lamb shift" . In addition, a correction as a result of distortion of 
the wave function of the 2s-state due to vacuum polarization, has recently been calculated 
by Ivanov et al. [68j. Numerically the results agree with previous calculations O [13] for 
the case of the hyperfine structure of the 2s-state (as given in Eq.[T9]), when the effect of 
finite nuclear size is neglected. 

Pachucki [1] has estimated an additional contribution of -0.005 meV corresponding 
to a vacuum polarization insert in the external photon for hydrogen. It is not clear 
to what extent this contribution is described by corrections to (V'^V) as described in 
Appendix C. A more recent calculation by Jentschura [57j gives a different result for this 
contribution. An independent recalculation is desirable and this will be presented for a 
part of the contribution in Appendix C. The effect of finite nuclear size (analogous to the 
Bohr-Weisskopf effect known for the hyperfine structure) will also be estimated there. 

The numerical value of the partial correction to the 2s-2p transition calculated there is 
of the order of 0.002 meV for hydrogen, which is too small to account for the discrepancy 
in the proton radius. It is not included in Table[51 since corrections to the Bethe sum 
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have not been calculated. However, this introduces an uncertainty in the theoretical value. 
Also the corrections are larger for other light nuclei. The total corrections calculated in 
Appendix C should be considered as uncertainties in the theoretical value for the Lamb 
shift. 

Summary of contributions 

Using the fundamental constants from the CODATA 2006 compilation ([E]), except 
where noted, one finds the transition energies in meV in Table [31 Here the main vacuum 
polarization contributions are given for a point nucleus, using the Dirac equation with 
reduced mass. Relativistic recoil corrections are given separately. Some uncertainties 
have been increased from the values given by the authors, as discussed in the text. 



Contribution 


Value (meV) 


Uncertainty (meV) 


Uehling 


zOo.Ozoz 




Kallen-Sabry 


l.OUol 




V JT ILeiaLlOllb lOOJ 


n 1 ^07 




sixth order [53] 


0.00752 




Total" LBL" [27] 


-0.00089 


0.00002 


mixed mu-e VP 


0.00007 




hadronic VP 


0.011 


0.001 


recoil [3] (eql36) 
recoil, higher order [3] 


-0.04497 
-0.0100 




recoil, finite size [H] 


0.013 


0.001 


recoil correction to VP [T] 


-0.0041 




additional recoil [HE] 


0.0575 




muon Lamb shift 






second order 


-0.66788 




higher orders 


-0.00171 




nuclear size (i?p=0.875 fm) 
main correction B ■ (r^) 




0.007 fm 


-4.002 


0.064 


Zemach moment [41j 
remaining order (aZ)^ [¥T] 


0.0244 


0.002 


-0.0001 




polarization 


0.0127 


0.003 


correction to the 2pi/2 level 


0.00004 





Table 3: Contributions to the muonic hydrogen Lamb shift. The proton radius is taken 
from |15j . 

In the case of the muon Lamb shift, the numbers in Table[3] are for the 2si/2 — 
transition. 

As calculated in appendix B (see Tabled?]), the finite size corrections for hydrogen can be 
parametrized as —5. 22718 (r^) + 0.00913 (r'^) (2), where energies are in meV and radii in fm. 
The total coefficient of (r^) differs slightly from that given in ref. [2] (-5.2262meVfm~^). 
The difference is due partly to a more precise determination of the coefficient of order 
(aZ)^ and partly to the inclusion of the radiative correction given by Eides and Grotch 
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Also, the shift of the 2pi/2-state is included. The second term can be approximated 
by 0.0365(18)((r2))3/2. This results in a total transition energy (without hyperfine inter- 
action) of 206.0592(60) - 5.22718(r2) + 0.0365(18) (r2)3/2. if the proton radius is taken 
to be 0.842(1) fm, the (total) nuclear size correction becomes —3.6855 ± 0.010 meV. With 
a radius of 0.875(7) fm, it is -3.978 ± 0.065 meV. 

A recent paper by Carroll et al. [SH] calculates the Lamb shift in muonic hydrogen rela- 
tivistically and nonperturbatively. This is very useful progress in removing the limitations 
of perturbation theory. For what is supposed to be the same calculated transition energy 
they obtain 206.0604 - 5.2794(r2) + 0.0546(r2)3/2. it would have been helpful if they 
had provided more details about how they obtain corrections due to two- and three-loop 
vacuum polarization contributions (Kaellen-Sabry, sixth order, and higher orders), muon 
selfenergy and muon vacuum polarization, hadron vacuum polarization, relativistic recoil, 
and nuclear polarization. Instead they simply used the values listed in the supplement to 
ref. |2]. Also, there are some problems with their calculation of the hyperfine splitting of 
the 2s state. 

Summary of contributions for muonic deuterium 

Muonic deuterium is in many ways similar to muonic hydrogen, but there are some 
differences. In addition to the different mass, the deuteron has spin 1 and both magnetic 
and quadrupole moments. For deuterium, one finds the transition energies in meV in 
TablelH Also here the main vacuum polarization contributions are given for a point 
nucleus, using the Dirac equation with reduced mass. 

As calculated in appendix B (see Table d^ . the finite size corrections for the 2s-2p 
transition in muonic deuterium can be parametrized as — 6.10940(r^) + 0.0112(r^)(2). 
where energies are in meV and radii in fm. The second term can be approximated 
by 0.0448(22)((r2))3/2. For more details, including the contribution for remaining or- 
der {aZy, see appendix B. The total transition energy (without hyperfine structure) 
calculated here is then 

230.2972 ± 0.04 - 6.10940(r2) + 0.0U8{22){{r^)f/^meV 

A very recent calculation of the same transition energy by Krutov and Martynenko [60] 
gives similar results, with one significant difference. For the nuclear structure contribu- 
tions of order (Za)^ they use the value given by Pachucki [52], which supposedly includes 
nuclear polarization and the Zemach moment term. In this work, the more standard sepa- 
ration is used. The total calculated here is approximately 0.433(21) + 1.50(3)=1. 933(45) meV, 
to be compared with 1.680(16) meV. The difference amounts to 0.256(60) meV, which is 
not negligible. An independent calculation of these nuclear structure effects would be 
very desirable. 
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-0.002001 




nuclear size (i?^^=2.130 fm |18j) 




0.003 fm 


main correction B • (r^) 


-27.718 


0.078 


Zemach moment [H] 


0.4329 


0.021 


remaining order (aZ)^ [H] 


0.0034 




polarization |50] 


1.50 


0.03 


correction to the 2pi/2 level 


0.00038 





Table 4: Contributions to the muonic deuterium Lamb shift. The deuteron radius is 
taken from [18j . 

Summary of contributions for muonic helium 

For ^He one finds the transition energies in meV in table |5l Also here the main vacuum 
polarization contributions are given for a point nucleus, using the Dirac equation with 
reduced mass. In the case of '^He results for two different radii (from [T^ I2U] ) are given. 
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0.0149 




polarization [H [53] 




3.1 


0.6 



Table 5: Contributions to the muonic helium Lamb shift. Finite size contributions are 
given for two values of the nuclear radius of ^He. 

As given in TablelT^in Appendix B, the finite size corrections for ^He can be parametrized 
as — 106.340(r^) -|- 0.400(r^)(2). where energies are in meV and radii in fm. The second 
term can be approximated by 1.40(4)((r^))^/^. The total transition energy would then be 
1667.937±0.005+3.1±0.6-300.477±1.4 +6.650±0.19 = 1377.210±1.5 meV. 

In previous work [H [8], recoil corrections to the Lamb shift in ^He denoted by "two 
photon" and "Breit" were given as — 0.44meV and -1-0.28 meV, respectively. These clearly 
correspond to corrections given here as (recoil p] (eql36)) of — 0.433 meV and (recoil, finite 
size IH]) of -0.1221 fm~^(r) (2) = 0.266 meV, respectively 

For ^He one finds the transition energies in meV in Table El Also here the main 
vacuum polarization contributions are given for a point nucleus, using the Dirac equation 
with reduced mass. From TablelT^ in Appendix B, the finite size corrections for ^He 
can be parametrized as —103. 508(5) (r^) -|- 0.3860(r^)(2). where energies are in meV and 
radii in fm. The second term can be approximated by 1.35(4)((r^))^/^. This gives 
a total value of 1648.402±1.3-400.075±3.67+0.015-f 10.258±0.305 = 1258.600±1.5 meV. 
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main correction B • (r^) 


-400.075 


3.67 


Zemach moment [H] 


10.258 


0.305 


remaining order [H] 


0.120 




correction to the 2pi/2 level 


0.0213 




polarization [T] 


4.9 


1.0 



Table 6: Contributions to the muonic helium Lamb shift. The nuclear radius of '^He is 
taken to be 1.966(10) fm. 

Here the estimated uncertainty does not include the uncertainty in the radius, but only 
the uncertainty in the coefficient of the Zemach moment and the nuclear polarization. 
This agrees quite well with a previous calculation (TU], which was considerably less precise, 
when the change in the contribution due the more recent measured radius is taken into 
account. 

Fine structure of the 2p state 

The fine structure of the 2p states can be calculated by using the relativistic Dirac 
energies, computing the vacuum polarization contributions with Dirac wave functions, 
and including the effect of the anomalous magnetic moment in the muon Lamb shift. An 
additional recoil correction (Eq JTU]) also has to be included. The results are summarized 
in TableE One should also include the i?^/2M7v-type correction to the fine structure, 
(see |3], Eq(38)). This is tiny (5.7 ■ 10~^meV for hydrogen) and is not included in the 
table. As mentioned before. Friar [H] has given expressions for the energy shifts of the 
2p-states due to finite nuclear size. A correction proportional to (r^) affects only the 2pi/2 
state, and thus contributes to the fine structure. The contribution to the fine structure 
of the 2p-state of hydrogen (4 ■ 10~^ meV) is negligible. However, this contribution is not 
negligible for the helium isotopes. For hydrogen, this result for the fine structure was 
subsequently reproduced by Martynenko [61]. His results for the hyperfine structure of 
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Frvdroffen 


deuterium 


^He 


^He 


Dirac 


8.41564 


8 86430 


144.416 


145.718 


IJehlinffl'VP') 




00575 


270 


275 


Kallen-Sabry 


0.00004 


0.00005 


0.002 


0.002 


anomalous moment 










second order 


0.01757 


0.01949 


0.3232 


0.3290 


higher orders 


0.00007 


0.00007 


0.0012 


0.0013 


Recoil (Eq.(IIO])) 


-0.08621 


-0.02521 


-0.1898 


-0.1107 


Finite size 


-0.00004 


-0.00027 


-0.0158 


-0.0119 


Total Fine Structure 


8.3521 


8.86419 


144.809 


146.203 



Table 7: Contributions to the fine structure {E(2p^i2) — E{2pi/2)) of the 2p-state in 
muonic hydrogen, deuterium, ^He and ^He, in meV. 



the 2p levels in hydrogen also agree with the results given below (and previously [H]). For 
deuterium, the fine structure agrees with that given by Krutov and Martynenko 



Hyperfine Interactions 

The Breit equation [3l O |58] contributions to the fine- and hyperfine interactions for 
general potentials and arbitrary spins were given by Metzner and Pilkuhn [62] • Here 
a version applicable to the case of muonic atoms {Zi = —1, si = 1/2, mi = m 



^l1 



"At; 



Z) is given. For most of the following, the potential is approximated 



Coul 



-aZ/r. 



V, 



L,si 



2m, 



IdV 
r dr 



Sirrir 



1 

m,. 



L ■ si 



(11) 



This can be rearranged to give the well-known form for spin 1/2 nuclei with an anomalous 
magnetic moment, namely 



1 dV 1 + + (a^ + l/2)mjq /m^ 
r dr iTii^m^ 



L ■ a„ 



Note that 



mNTTlf, 



+ 



2ml 



2ml 



2m% 



so that the terms not involving in the spin-orbit contribution are really the Dirac fine 
structure plus the Barker-Glover correction (Eq. [lU]) . 
Also 



1 1 dV 

2/712 ^ 



1 + ^2 



S2mr 



1 

m2 



L ■ S2 



Usually one writes 



Z(l + K2) 



mr. 



777-2 iii'p 

where H2 is the magnetic moment of the nucleus in units of nuclear magnetons [fi^ = e/2mp) 
A value of fid = 0.85744 fj,^ = 0.307012 /ip corresponds to Kd = 0.714. In the case of ^He, 
the magnetic moment is fiHes = —2.127b resulting in nuez = —4.184 
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The spin-spin interaction is given by 

_ 2(l + a^)/i2 rl dV 



2s2m^m2 



r dr 



(3Si ■ fS2 ■ f - Si ■ S2) - -V^V^Si ■ S2 



For deuterium the quadrupole moment also contributes, with 

Q 1 dV 



Vn 



2 r dr 



3S2 ■ fS2 ■ f - S2 ■ S2 



with Q in units of l/m|. The quadrupole moment of the deuteron is taken to be 
Q = 0.2860(15) fm^ [ISlESlEZ]. In other units, one finds Q = 25M/ml = 7.345x10-^ MeV'^. 

The Uehling potential has to be included in the potential V{r). For states with i > 
in light atoms, and neglecting the effect of finite nuclear size, we may take 



IdV 
r dr 



aZ 



1 + 



2a 
3^ 



00 .2 



1) 



1/2 



1 

2^2 



'1 + 2merz) ■ e 



-2meTZ 



dz 



(12) 



which is obtained from the Uehling potential [l]|29l|30] by differentiation. Then, assuming 
that it is sufficient to use nonrelativistic point Coulomb wave functions for the 2p state, 
one finds 

■ldV\ /ldV\ 
^r dr / 2p \r dr / 2p 

where 



2a 
37r 



00 / 2 



1) 



1/2 



1 + az) 



+ 



2az 



az] 



dz 



(13) 



with a = 2me/{aZmr). For hydrogen, £:2p = 0.000365, for deuterium £2^ = 0.000391 and 
for ^He e2p = 0.000894. 

The correction due to vacuum polarization (Eq. (|T3|) ) should be applied to the HFS 
shifts of the 2p-states, (and also to all spin-orbit terms in the "muon Lamb shift"). 

Note that for the 2p state 



ldV\ 
r dr I 



2aZ- 



{aZrrir/nY 
^+1)(2^+1) 



aZ{aZmr 
24 



The hyperfine structure can also be calculated relativistically, using the formalism 
given in [1]. The finite extent of the magnetization density and the effect of vacuum 
polarization should be taken into account also in this approach. 



Hyperfine structure of the 2s-state: 

The expectation value of ^^^^2 in an n-s state with j = 1/2 is 

Ai?„, = = ^ff"^^'""' ■ (1 + a,) [F{F + 1) - S2{S2 + 1) - 3/4] 
3n-^m^m2S2 

When S2 = 1/2, and ^2/^p = (1 + '^2)/^2i this reproduces the well-known result for 
muonic hydrogen (see, for example [3], Eq. (271,277): 

AE„, = -^.(l + K2)-(l + a^)(VV)[F(F + l)-3/2] = ^ . (1 + a^)[F(F + 1) - 3/2] 
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with ^ ^ 

/3 = ■ (1 + t^2) = i8/n^) X 22.8054 mel^ (14) 

The numerical value was calculated for hydrogen. Note that /3 ■ (1 + a^) = 22.8320 meV, 
which is the well-known value. Since ^He also has spin 1/2, the same formula is valid, 
and one obtains /3 = -171.3964 meV. (/3 ■ (1 + a^,) =-171.5963 meV) 
For deuterium, with S2 = 1, the corresponding hyperfine splitting is 

AE^, = ^MM . (1 + . (1 + a,) ■ [FiF + 1) - 11/4] 

= ^ ■ (1 + a^) ■ [F{F + 1) - 11/4] = ■ 2.04766 meF ■ [5^,3/2 - 2Sf,i/2] 

for a total splitting of 6.14298meV in muonic deuterium. This is in reasonably good 
agreement with the result given by Carboni [T2] . 

The QED corrections to the energy shift of the 2s-state were discussed in [3l [9] . In 
muonic hydrogen they are given by: 

AE2s = /3 ■ {l + ttf,) ■ {l+evP + e^ertex+£Breit + £Zem) ■[Sf1-3Sfo]/4: (15) 

The corrections due to QED effects, nuclear size and recoil are analogous for muonic 
deuterium and ^He. Other corrections due to recoil, nuclear polarization, the weak 
interaction and so on can be included by adding corresponding e's. Here ([3], Eq. (277)) 

17(aZ)2 _4 

^Breit — 3 — 1.13-10 ■ Z 

o 

is a relativistic correction that gives the difference between the value obtained in rela- 
tivistic perturbation theory using Dirac wave functions, and the nonrelativistic value (to 
leading order in (aZ)^). The vertex correction (|3l |65]) is given by 

e.ertex = a{Z o) {\n{2) - = -0.9622 • 10"^ • Z 

This includes a correction of 3a(Za)/4 = 0.3994 ■ 10~^ ■ Z due to muon loop vacuum 
polarization. Corrections due to hadronic vacuum polarization can be expected to be 
comparable to this value. Some higher order corrections to Svertex of order a{ZaY \v?{Za) 
are possibly numerically important and were given in Ref. |65j as 



(-^ ln'((^«)"') + C22(2) ln((Za)-2) J = -0.0073 ■ 10"^ {Z = 1) 

This correction adds an additional -0.00017 meV to the hyperfine splitting in the case of 
hydrogen (and -0.00004 meV in muonic deuterium). If Z=2, the numerical value of the 
extra correction is —0.211 x 10~^ and the contribution to the hyperfine splitting in muonic 
^He is 0.0036 meV. 
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The main vacuum polarization correction has two contributions. One of these is a 
result of a modification of the magnetic interaction between the muon and the nucleus 
and is given by (see [10] ) 



£vpi 



4a r fRrJr) 
r ar 



37r2 



Rns{0) 



CO ^^2 _ ^)l/2 



1 + ^ 



(fjQ{(ir)GM{q) dq 



dz 



(16) 



2zy Aml[z'^ + {q/2meY] 



One can do two of the integrals analytically and obtains for the 2s-state (with a = 2me/ [aZrUr 
and sinh(0) = g/(2me) = K/a) 



£vpi 



4a 
3^ 



F{(P)GM{aZmrK) dK 



7 



+ 



6 



(1 + ^2) (1 + 7^2)2^ 



(17) 



where F{(j)) is known from the Fourier transform of the Uehling potential and is given by 
EqdH). 

The other contribution, as discussed by [651 [66] arises from the fact that the lower 
energy hyperfine state, being more tightly bound, has a higher probability of being in a 
region where vacuum polarization is large. This results in an additional energy shift of 

2 / Vuehi{r)ip2s{r)SMip2s{r)d\ 



13?/ - 3 - 2y2 1 ■ 

4 Ay 
(18)- 



Following Ref. [65] with y = {aZmr/2) ■ r, one has 



aZnir 



(l-y)(ln(2y) + 7) + 



(7 is Euler's constant), and 

^2s(r) = V2s(0)(l-2/)e-^ 
One finds after a lengthy integration 

evP2 = ^ r T^^oGE{o^ZmrK)F{<P) 



37r2 7o 1 + 7^2 
1 



17 



+ 



41 



24 



+ 



2 (l + i^2)2 (l + i^2)3 (1 + 7^2)4 

ln(l + 



1 + ^2 



7 



6 



+ 



K 



1 



(l + ir2) (1 + 7^2)2 

19 20 



+ 



12 



2(1 + 7^2) (1 + 7^2)2 (1 + 7^2); 



(19) 



Sternheim[5B] denotes the two contributions by 5m and 5e, respectively. An alternative 
exression, obtained by assuming a point nucleus, using Eq.(131) from [Ij for the Uehling 
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potential, and doing the integrations in a different order, is 



£VP2 — 




az I 23 3 



(20) 



2 1 + az 8(l + az)2 2{1 + azf 



+ ln(l + az) ■ ( 1 - — + — ^ — 1 

^ ^ V l + az 2{l + azyj 



dz 



with a = 2me/ (aZmred) as in the case of e2p- Both methods give the same result when 
the effect of nuclear size is neglected. In the case of ordinary hydrogen, each of these 
contributes 30:^/8 = 1.997 ■ 10~^. The accuracy of the numerical integration was checked 
by reproducing these results. One can thus expect that muonic vacuum polarization will 
contribute 3a^/4 ~ 4 ■ 10~^, as in the case of normal hydrogen. This correction was 
included in e^ertex- Martynenko [70] includes only one of these contributions. The energy 
shift is approximately 0.00091 meV in hydrogen. Contributions due hadronic vacuum 
polarization should be slightly smaller. An estimate following the prescription given in 
[39] gives a correction of approximately 0.2667 ■ 10~^ ■ Z, or approximately 0.0006(1) meV 
in hydrogen. Contributions due to the weak interaction were estimated by Eides [75] 
and give 5.8 x 10~'^(m^/me)-EF which works out to 0.00027 meV for the weak interaction 
contribution in muonic hydrogen. Unfortunately this calculation does not seem to be 
directly applicable to other nuclei. 

Finite nuclear size has an effect on the value of eyp- For muonic hydrogen, one obtains 
£ypi=0. 00211 and eyp2=0. 00326 for a point nucleus and ev/pi=0. 00206 and ei/P2=0. 00321 
with a proton (magnetic) radius of 0.875 fm. Including this effect reduces the total 
contribution due to vacuum polarization by 0.00023 meV. 

For muonic ^He, one obtains eypi=0. 00295 and ev/P2=0. 00486 including the effect of fi- 
nite nuclear size. For a point nucleus the values would be £:ypi=0. 00315 and £v'P2=0. 00506. 

For the case of muonic deuterium, £ypi=0. 00218 and £v/P2=0. 00337 for a point nucleus. 
Including the effect of finite nuclear size gives £i/pi=0. 00207 and £yp2=0. 00326. 

The contribution to the hyperfine splitting of the 2s-state of muonic hydrogen (with 
finite extension of the magnetization density) is 0.04703 me V+0. 07328 me V=0. 12031 meV. 
The combined Breit and vertex corrections change this value to 0.11794 meV. 

The contribution to the hyperfine structure from the two loop diagrams [32] can be 
calculated by replacing U2{aZmrK) = {a/3'7i)F{(j)) by U^^aZnirK) (as given in [U |7]) 
in equations [17] and [191 The resulting contributions in hydrogen are 1.636 ■ 10~^ and 
2.460 ■ 10~^, respectively, giving a total shift of 0.00093 meV. Martynenko [70] neglected 
the contribution corresponding to Eq.[Tni The contributions in muonic deuterium are 
1.688 ■ 10~^ and 2.545 ■ 10~^, respectively. For muonic ^He they are 2.511 ■ 10~^ and 
3.928 ■ 10-^ respectively. 

The correction due to evP2 can be regarded as a result of distortion of the wave function 
of the 2s-state, and has recently been calculated by Ivanov et al. pS]. Numerically their 
results agree with those given here to within a few one per cent, for both the second order 
and fourth order contributions in hydrogen and in deuterium. A more complete calculation 
of the vacuum polarization corrections to the 2s-state in hydrogen by the same authors 
[69] gives the same result for all these corrections to three significant figures, when nuclear 
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size is neglected. This paper includes the contribution of three further graphs, which give 
an additional contribution to eyptot of 1.64 ■ 10~^, and a total two-loop contribution of 
0.00130 meV to the hyperfine splitting of the 2s-state of hydrogen. 

The main correction due to finite extension of the magnetization density is known as 
the Bohr-Weisskopf effect (^67^) and is equal to 

ezem = -2aZm^(r)(2) 

where {r)(2) is given in p, I^Tj H5]. 

Ref. |1] claims that this correction does not treat off-shell effects and/or recoil properly. 
For hydrogen, Carlson et al. [73] give a recoil correction corresponding to e^ec = 0.00093, 
resulting in a correction of 0.02123 meV. 

For hydrogen the value of 2aZmr is 0.007024/m~^. Using the value 
(^)(2) = 1-086 ± 0.012 fm from [43j, gives a Zemach correction of ezem = —0.00762, and 
a contribution of -0.1742(19) meV to the hyperfine splitting of the 2s state in hydrogen. 
Distler et al. [H] give {r)(2) = 1-045 ± 0.004 fm for the 'magnetic' Zemach radius; the 
resulting value of ezem is -0.00734, and an energy shift of -0.1676(6) meV. 

The corrections due to finite size and recoil have been given in [1] as -0.145 meV, while 
a value of -0.152 meV is given in |70]. Combining the recoil correction from Ref. [73] with 
the Zemach correction results in values compatible with either of these, depending on the 
Zemach radius (for example, -0.1464(6) meV when the Zemach radius of Distler et al. [H] 
is used). 

A correction for possible nuclear polarization effects has been calculated by Chered- 
nikova et al. [71] with the result Epoi = 0.00046(8), for an additional contribution 
of 0.0105 ± 0.0018 meV. Carlson et al. [73] give a value of Spoi = 0.000351(114), or 
0.0080 ± 0.0026 meV for this correction. To be consistent with the recoil correction given 
above, the correction of Carlson et al. should be used (see also ref. [71]). 

Comparable calculations for muonic deuterium and ^He have not been found. 

In hydrogen, the total corrections (other than for the Zemach radius) amount to 
0.1785 ± 0.0027meV (or 0.1762 ± 0.0027meV if the finite extent of the magnetization 
density in the VP correction is included). Adding Ep gives 22.9839±0.0027meV (or 
22.9816±0.0027meV). The correction due to the Zemach radius is -0.16037meVfm~^(r)(2). 
If the value of the Zemach radius from Distler et al. jH] is used, the total hyperfine 
splitting in the 2s-state of muonic hydrogen is 22.8140 meV. 

Alternatively, one could replace the contribution from the Zemach correction by a 
polynomial expansion in the magnetic radius, as suggested by Carroll et al.[59]. Since 
the leading term in a perturbative calculation (relativistic or nonrelativistic) is linear, 
this should take the form (ar -|- 6r^). Since these authors take most of their radiative 
corrections from other work [TOj, which neglects some of the corrections included here, a 
more detailed comparison is not appropriate. 

For muonic deuterium. The coefficient of (^')(2) is -0.007398 fm~^, giving, with 
(r)(2) = 2.593 ± 0.016 fm from p], ezem = -0.01918 ± 0.00012. Nuclear polarization 
and recoil corrections may be important, but have not been calculated. The total hyperfine 
splitting of the 2s-state, including all known corrections, is 

3 

AE2s = -/3d ■ (1 + fl/x) ■ (1 + £VP + ^vertex + ^Breit + £Zem) = 6.0584(7) UieV 
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For muonic ^He, The coefficient of (?")(2) is -0.01506 fm~^, giving, with (r)(2) = 2.403 fm 
(a Gaussian charge distribution was assumed), Ezem = —0.0362. The total hyperfine 
sphtting of the 2s-state, including all known corrections, is -165.7949 me V. 



Hyperfine structure of the 2p state 

The hyperfine structure of muonic hydrogen is calculated in the same way as was done in 
earlier work [HI HO] , but with improved accuracy. Most of the formalism and results are 
similar to those given by [1] and [63] . 



Hydrogen 

The hyperfine structure of the 2p-state in hydrogen is given by [HI [63] {F is the total 
angular momentum of the state) 



1 



/ldV\ 



4:771 dr ' 2p 

+ Q33'{CFi{l + a^)-2{l + x)) 



2(1 + x)5,,,[F{F + 1) - j(j + 1) + 3/4] 



where j = ^2j + 1, the 6-j symbols are defined in 
Cpi = 5fi - 25fo - (1/5)5f2 
Also 



F ] 

1 . 

2 ^ 

and 



F 
1 

2 
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X 



m^(l + 2k) 



2mAr(l + K) 

represents a recoil correction due to Thomas precession [9l[58l[63j. For muonic hydrogen, 
X = 0.0924. 

As has been known for a long time [U M, [101 [M], the states with total angular 
momentum F = 1 are a superposition of the states with j = 1/2 and j = 3/2. Let 
the fine structure splitting be denoted by 5 = E2pz/2 — -^2^1/2, and let (3' = (3 ■ {1 + e2p), 
to take the correction due to vacuum polarization into account. (/3 was defined above, in 
Eq.M)- For hydrogen, e2p = 0.000365 and /?' = 22.8138 meV 

The energy shifts of the 2p-states with total angular momentum F (notation '^^^^Lj) 
are then given in Table [S] 



State 


Energy 


Energy in meV 


^Pl/2 


-P'{2 + x + a^)/8 


-5.9704 


^Vl/2 


(A - R)/2 


1.8458 


'^P3/2 


(A + R)/2 


6.3760 


^P3/2 


6 + + 5x/4 - a^/4)/20 


9.6242 



Table 8: Hyperfine structure of the 2p-state in muonic hydrogen. Here 5 
is the fine structure splitting of the 2p state. 



^.352 meV 



where 

A = 6-/3'{x-a^)/16 
= [6- /3'(1 + 7x/8 + a^/8)/6]' + {f3')\l + 2x - a^f/2%% 
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Some minor errors in [9j have been corrected. These numbers differ shghtly from those 
given in ref. |3]. 

Helium-3 

The formulas for muonic ^He are identical to those for hydrogen, but the numerical values 
are, of course, different. For muonic ^He, x = 0.0435, 6 = 144.809 meV, /3 = -171.396 meV 



State 


Energy 


Energy in meV 




-/3'(2 + X + a^)/8 


43.8458 




(A - R)/2 


-14.7877 




(A + i?)/2 


160.0510 


^P3/2 


6 + /3'{1 + 5x/4 - a^/4)/20 


135.7673 



Table 9: Hyperfine structure of the 2p-state in muonic '^He. Here 6 = 144.809 meV. 

£2^=0.000894, and = /3 ■ (1 + £2?) = -171.550 meV. The energy shifts of the 2p-states 
with total angular momentum F (notation "^^^^Lj) are then given in TableO 

Deuterium 

For the 2p state, the matrix elements of the magnetic hyperfine structure have been given 
by Brodsky and Parsons [63]. For hydrogen they are the same as those calculated above. 
The matrix elements for the magnetic hyperfine structure are then given by 



j 


3' 


Energy 


1/2 


1/2 


{I3d/Q){2 + xa + a^)[-6F,i/2 + 1/25^,3/2] 


3/2 


3/2 


6 + (/3d/4) (4 + 5xd-a^) [-1/6 5^,1/2- 1/15 5i.,3/2 + 1/10 5f,5/2] 


3/2 


1/2 


(/3z5/48)(l + 2xd - a^)[-V26p,i/2 - V56f,3/2] 



where was defined in the previous section. Numerically it is equal to 4.0906 meV. 
Xd = {ml/mdmr)iKd/il + Kd)) = 0.02A8. 

For the evaluation of the contributions of the quadrupole HFS, let 

^ g /I 

2\r dr / 

Numerically one finds, for a point Coulomb potential, and the 2p-state, 

eo = = 0.43423 meV^, 

^ 2 24 

As mentioned before, the Uehling potential has to be included in the potential V{r). 
For states with £ > in light atoms, this can be taken into account by multiplying 
and eg by (1+ £2p) where e2p is given by Eq. (fT3|) . With a numerical value of e2p = 0.000391 
for muonic deuterium, the value of eg is increased to 0.43439 meV and the value of is 
increased to /3^=4.0922meV. 
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The quadrupole interaction results in energy shifts of 



j 


f 


Energy 


1/2 


1/2 





3/2 


3/2 


eg [Sf,i/2 - 4/5 5f,3/2 + 1/5 5f,5/2] 


3/2 


1/2 


eg [y/2 5f,1/2 - l/VS 5f,3/2] 



Then for the 2p-level with j = j' = 3/2 and F = 5/2, the energy shift is given by 
5 + eg/b + (/3d/40) (4 + 5xrf-a^) = 9.3728 meV. For the 2p-levels with F = 1/2 and 
F = 3/2, the corresponding matrices have to be diagonahzed. The resulting numerical 
values for the eigenvalues are given in TablefTUl 



State 


Energy in meV 


Vl/2 


-1.4056 




8.6194 




0.6703 


^P3/2 


8.2560 


^P3/2 


9.3728 



Table 10: Hyperfine structure of the 2p-state in muonic deuterium. Here 
5 = 8.86419 meV. 

Table [TT] gives the contributions to the transition energies due to fine and hyperfine 
structure in muonic hydrogen relative to the 2s-2pi/2 transition energy given in TableEi 



Transition 


Energy shift in meV 


Vl/2 Si/2 


-11.6739 


^Pl/2 Si/2 


18.9563 


V/2 Si/2 


-3.8577 


^P3/2 Si/2 


23.4865 


^P3/2 Si/2 


0.6725 


^P3/2 Si/2 


3.9207 



Table 11: Fine- and hyperfine contributions to the Lamb shift in muonic hydrogen. 

Table [T^] gives the contributions to the transition energies due to fine and hyperfine 
structure in deuterium. 

TabledS] gives the contributions to the transition energies due to fine and hyperfine 
structure in muonic ^He, relative to the 2s-2pi/2 transition energy given in TableEl 
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Transition 


Energy shift in meV 


I-'l/z '^1/2 


2 6333 


P3/2 ■'1/2 


12 6583 


4 2 


4.7092 


V3/2 Si/2 


12.2949 


^Pl/2 Si/2 


-3.4251 


^P3/2 Sl/2 


6.5999 


4 4 
>l/2 - Si/2 


-1.3492 


4 4 
P3/2 - Si/2 


6.2365 


6 4 
T3/2 - Si/2 


7.3533 



Table 12: Fine- and hyperfine contributions to tlie Lamb shift in muonic deuterium. 



Transition 


Energy shift in meV 


Vl/2 Si/2 


85.294 


Si/2 


-139.134 


^Pl/2 Si/2 


26.661 


^P3/2 Si/2 


35.705 


^P3/2 Si/2 


201.500 


^P3/2 Si/2 


177.216 



Table 13: Fine- and liyperfine contributions to the Lamb shift in muonic ^He. 



Conclusions 

Experimental precision has reached the point that some previously neglected effects should 
be taken into account. A few of these have been discussed in this paper 

The consequence of the model dependence of the coefficient of (in units meVfm"^) 
for the determination of the nuclear radius has been discussed. The conversion of the 
coefficient of the third Zemach radius to ((r^))^/^ is not unique. For the case of hydrogen, 
the contribution 0.009126 meVfm"^ ('^^)(2) is most plausibly 0.0365(18) meVfm~^((r^))^/^. 
This will increase somewhat the theoretical error given in [2J. Assuming the value 
of the proton radius given in [2], the contribution to the transition energy would be 
0. 02179±0. 00107 meV (here only the uncertainty in the coefficient was considered). 

The radiative correction to the contribution to the energy shift of the 2s-state due to 
nuclear size has been included and the nuclear size contribution to the vacuum polarization 
correction has been calculated more completely. Details are given in Appendix B. 

The previously neglected corrections to the muon Lamb shift discussed in Appendix C 
have not been included in the summaries. They introduce an additional theoretical 
uncertainty, which is of the order of 1-2 /leV for hydrogen. They are not negligible for 
the other cases studied here. Although the basic conclusions in the hydrogen experiment 
are unchanged, the theoretical uncertainty is increased, which would slightly increase the 
uncertainty in the determination of the proton radius. 

The numerical results for muonic helium differ somewhat from those of Martynenko 
j76] . It is difficult to compare all terms, but Martynenko's calculation of the Wichmann- 
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Kroll contribution is incorrect (the sign is wrong, for reasons described in section 9.3.1 
of ref. [3j). Also, not all of the fourth order contributions to the muon Lamb shift were 
taken into account. 

The hyperfine splitting of the 2s-state has been recalculated, including some previously 
neglected effects. The result for muonic hydrogen is 

(22.9816±0.0027)meV-0.16037meVfm-i(r)(2). 

For muonic deuterium the total hyperfine splitting of the 2s-state is 

6. 17625 meV-0.045446 meV fm^^ (r) (2) = 6.0584(7) meV 

A relativistic nonperturbative calculation of these effects is desirable, and has been 
performed in part [59]. Since some problems with the hyperfine structure of the 2s state 
have been mentioned, no further comments will be given here. 
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Appendix A: Details of the relativistic recoil calculation 

As mentioned before, the energy levels of muonic atoms are given, to leading order in 
l/rriN by 

where E^. is the energy level calculated using the Dirac equation with reduced mass and 
Bq is the unperturbed binding energy. Also 

h{T) = -Pi(r)(Pi(r) + -Q,ir)) - ^Q,(r)[P,ir) + ^] 

where Pi, (52, and Q4, are defined in Eq.([9]). The calculation of these corrections in first 
order perturbation theory, and to first order in the Uehling potential, is described in detail 
here. It was also given in ref. [6]. Some minor mistakes in the numerical calculation have 
been corrected. In the case of a point Coulomb potential {y{r) = —aZ/r), the functions 
Pi and Q4 are identically zero. The total potential considered here is 
V{r) = -aZ/r + Vuehiir). 

For the contribution due to the term 3^/2771^ the unperturbed binding energy is taken 
to be equal to Bqd + {Vuehi), where Bqd is the point Coulomb Dirac binding energy. The 
contribution to Bl/2mN that is linear in the Uehling potential is then {—BQj:,{Vuehi))/'mN- 

Note that 

Thus it is possible to rearrange h{r) to give 

h{r) = P,{r)V{r) - ^Q,{r)[P,{r) + 



3r r-^ 
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As noted elsewhere in this paper (see Eq.f ll6p ). an effective effective charge density 
pvp for vacuum polarization can be derived from the Fourier transform of the Uehling 
potential. Recall that (for a point nucleus) 

V„e.(r) = . X.(2m.r) 



Svr J I \ ) \TT Jq + Am'lz'^ 

where Xn{x) is defined in [1]. In momentum space, the Fourier transform of V'^V is 
obtained by multiplying the Fourier transform of V by —q^. One then obtains 

ATcpvp{r} = —■ dz 1 + 



Stt J I z"^ \ 2z^ ) \7r g2 + 4m^2;2 

?^f^2(g)jo(gr)rfg 



9 

^ ' ,2 



TT 



(22) 







U2{(l) is defined in p[] (see also Eq.H]). 

Keeping only the Coulomb and Uehling potentials, one finds 

2« 

Pi(r) = -aZ— (2me)xo(2mer) 

OTT 

Q2(^) = aZ ( 1 + — [xi(2mer) + (2mer)xo(2mer)] 



2\ 1"°° 1 (6gr — (gr)^) cos(gr) + (3(gr)2 — 6) sin(gr) 



) Jo + ^m'lz'^ q ^ 

Since vacuum polarization is assumed to be a relatively small correction to the Coulomb 
potential, it will be sufficient to approximate Q2{r)/r by aZ/r and 

Pi(r)(Pi(r) + -Q2{r)) = Pi{r)V{r) by -^Pi(r). 

Since the correction is to be calculated to linear order in the vacuum polarization 
potential, it will also be sufficient to use point Coulomb wave functions to calculate the 
expectation values. Also, for this case, one can use Schroedinger wave functions, since 
an accuracy of 1% will be adequate. After some algebra, one can reduce the expectation 
values to single integrals: 

(az)'^ — az + 1 ^ 1 



5m + 7^. T^Sn dz 



'1 + az)^ ™ (1 + az)^ 



, ^,3 2a r (^2. 1)1/2 / 1 
r 6t: Ji z \ 2z^ 

'2{azY + l 



N.^^oH ; Ti^n dz 

2{l + azY 2{l + azY 
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with a = 2me/{aZmr). 
Finally, 

^aZ^ , 2a i^^ - 1)'^' 



6 + 2, , 3{az)^ + 2az-l 
(- \n{l + az)) + 



+ az az 



+ 



3 + 

4(1 + az)4 



l-3az-2(a2)2 
+ 4(1 + az)4 



'X^azf 
dz 



Combining these expectation values according to equations [7] and |8|, one finds a contri- 
bution to the 2p-2s transition of -0.00414 meV (hydrogen) and -0.00246 meV (deuterium). 
To obtain the full relativistic and recoil corrections, one must add the difference between 
the expectation values of the Uehling potential calculated with relativistic and nonrel- 
ativistic wave functions, giving a total correction of 0.0169 meV for muonic hydrogen. 
This is in very good agreement with the correction of .0168 meV calculated by Veitia 
and Pachucki |56J. The treatment presented here has the advantage of treating the main 
contribution relativistically and applying a small correction that can be calculated using 
first order perturbation theory. For hydrogen, a very small correction of 0.0002 meV for 
the Kallen-Sabry contribution to the main term was included. For deuterium, one obtains 
a total correction of 0.0214 meV. 

A similar relativistic recoil correction for finite nuclear size should be included in a 
relativistic calculation starting from the Dirac equation (see Ref. [H]). However it is 
automatically taken into account when the calculation is based on a generalized Breit 
equation, such as in Ref. [56j . 



Appendix B: Higher order contributions to the correction for finite nuclear 
size 

If the transition energy is written in the form 



LS 



A + B{r^)+C{{r^)f/^ 



it is necessary to calculate the quantities A, S, and C . Suggested values for C have been 
given in the text, and A can be determined from the summary tables. Here the higher 
order contributions to B mentioned previously are explicitly calculated for all four cases. 
There are several contributions and one finds 

B = ha + h + hc + hd + he 
From Eq.O the main term in the expression for the energy shift of the 2s-state is 



2aZ ( aZiTLr 



n 



This defines ha- 
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Radiative corrections to the main term have been calculated by Eides and Grotch 
They contribute an additional correction of order ^(aZ)^, which amounts to 



laZ ( aZnir 



n 



a2Z(23/4-41n(2) -3/4)(r2) = ^^^(S - 4 ln(2)) ■ ^^(r^) = h{r^) 



Next, more details for the evaluation of the order {aZY contributions to the correction 
for finite nuclear size to the energy of the 2s state are given. 

The terms given by Friar [H] involving {aZ)^{FjiEL + ^rFNR) in equation (|5]) are given 
for n=2 explicitly by 



35 1 1 



and 



Fnr = l{{r'mi-l+Hc^Z)] + ^{r'){r'Hmrr))+^^^ 

The expectation values for the various moments of the charge distribution are given in 
Appendix E of ref. [H]. Approximate values for the integrals I^^^, /-f^^, -^2^^ and I^^ 
are given only for uniform and exponential charge distributions, so these extra terms, with 
the exception of the contribution to Frel that are proportional to (r^) were calculated for 
an exponential charge distribution, even if this is not completely realistic. The additional 
contribution proportional to (r^) in the finite size correction to the 2s level in equation (|5]) 
is 



2{aZY 



( aZvtirX ,9,. 35 , , 



The remaining terms are small, but a model independent evaluation is prohibitively 
difficult. Numerical values are given in Tables (3,4,5,6) with the heading "remaining 
order {aZf'\ 

There are two corrections due to finite nuclear size to the vacuum polarization contri- 
bution. One of these was obtained as a result of numerical integration of the expectation 
value of the Uehling potential, as discussed previously. It will be denoted by hd{r'^) and 
numerical results are listed in Table[TU The (small) correction arising from the Kallen- 
Sabry potential was also included. The other correction, due to iterations, has been 
calculated for hydrogen by Eides et al. [3] and by Pachucki [1]. It is given by 

AE2S = — (r^) / Vuekiir)iJ2sir)G',,ir,0)^,Md'r 



Ivanov and Karshenboim ^77j have given an expression for G"2g(r, 0). With x = aZnirr = br, 
it is 

G'^^{r, 0) = ' [4x(2 - x)(ln(a;) + 7) + ISx"^ - 6x - - 4] 
(7 is Euler's constant). Recall that 

i^2s{r) = V^2s(0)(l-x/2)exp(-x/2) 
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Then 



V'a.WG'a.lr, 0)^2.(0) = aZml^{l-x/2)e 



(1 -a;/2)(ln(a;) +7) + 



13a; -6 



X 



1 " 

2x 



It is possible to show that the integral is proportional to the correction evp2 to the 
hyperfine splitting of the 2s-state. If one rewrites Eq. (fT5]) in terms of a; = 2y, one finds 
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■i-K 



ilJ2s{r)5M^2s{r) = — T {aZm'^VF){l - x/2)e 



x/2)(ln(x) +7) + 



13x — 6 



X 



8 



1 

2x 



Thus 



and hence 



^2s(r)G^,(r,0)7/;2.(0) 



2s 



AnaZ , 



1 



3 'SnAup 
By definition, the quantity evp2 is given by 

2 



87r Aup 
Vuehi{r)^2sir)SM^2sir)d^r. 



£VP2 



Aup 



Vuehi{r)^lj2sir)6M^2sir)d r. 



Hence 
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2s 



aZ 
'I2" 



Since the two-loop VP corrections to evp2 have been calculated, they are included in 
be- Effectively, the one-loop correction is multiplied by approximately 1.0077. TabledH 
gives numerical values of these coefficients (in meVfm"^). The total for hydrogen differs 
slightly from that given in ref. [2] (-5.2262meVfm~^). The difference is due partly to a 
more precise determination of the coefficient be and partly to the inclusion of the radiative 
correction (6^). 





Hydrogen 


deuterium 


3He 


^He 


ba 


-5.1973 


-6.0730 


-102.520 


-105.319 


K 


- 0.00062 


-0.00072 


-0.0243 


-0.0250 


he 


- 0.00181 


-0.00212 


-0.1275 


-0.1310 


bd 


- 0.0110 


-0.0130 


-0.3176 


-0.3297 


be 


- 0.0165 


-0.02062 


-0.5217 


-0.5392 


52,=total (2s) 


-5.22723 


-6.10946 


-103.511 


-106.344 




-0.0000519 


-0.0000606 


-0.00409 


-0.004206 



Table 14: Contributions to the coefficient of (r^) for the energy shift of the 2s and 
2p-states in muonic hydrogen, deuterium, ^He and ^He, in meVfm"^. 



It is interesting to note that the total coefficient of (r^) calculated here for *^He is 
in excellent agreement with the calculated value given in ref. where the finite size 
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correction for ^He was estimated to be -106.2 (r^) + 1.4 (r^) (in meV). An additional term 
proportional to the Zemach radius is discussed in the section on relativistic recoil. For 
^He it was parametrized as 0.4 (r). It is actually proportional to (r)(2), and was evaluated 
in the section on relativistic recoil. 

For ^He the binding energy of the 2pi/2-state is decreased by 0.0148 meV (for a radius 
of 1.676 fm and a Gaussian charge distribution, which is a fairly good approximation for 
helium, according to Friar The contribution to the fine structure is 0.0118 meV. The 
difference is due to the term proportional to (r^), which is the same for both 2p-levels. 
If the radius is 1.681 fm both of these numbers are increased by 0.0001 meV. For ^He the 
binding energy of the 2pi/2-level is decreased by 0.0213 meV (for a radius of 1.966 fm). 
The contribution to the fine structure is 0.0158 meV. 

For muonic deuterium, the first term of equation (E]) would contribute 
-6.0730 (r^) =-(27. 817±0. 078) meV (using the value of the radius from the newest CO- 
DATA compilation [H]). If the radius of 2.130(3) from ref. [18] is used, this contribution 
is -26.8585±.076meV. Using the total coefficient from Table[Tl gives -27.718±.076meV. 

The term proportional to {r^)(2) (the coefficient is 0.0112 meVfm"^) gives a contri- 
bution of 0.382 meV or 0.417 meV, depending on the model for the charge density, and 
a radius of 2.14 fm. In ref. [13] it is suggested to calculate this term according to the 
prescription 

(r3)(2) ^ 4.0(0.2)((r2))^/2 

Using this value, one obtains 0.439(22) meV with a radius of 2.14 fm (0.433(22) meV with 
a radius of 2.13 fm). The last terms in Eq.(|5]) contribute a term proportional to (r^) and 
remaining terms. The first is bc{r'^) = -0.00966 meV; and the remaining terms (of order 
(aZ)^) given in [41] contribute 0.00337 meV. The complete contributions are given in 
TableU 

Appendix C: Further corrections 

A number of the corrections described here, in particular the "muon Lamb shift" and the 
hyperfine structure of s-states, involve the expectation value of V'^V. Note that using the 
normalizations of [H |9], one has V'^V = —AnaZp where p is the nuclear charge density. 

Usually V'^V = —AiraZ p{r) is approximated by a delta function, giving V = —aZ/r 
and 

(VV) = -A7caZ\4jUO)\%o = -AZa{^^f6,o 

However, the potential should be corrected for vacuum polarization due to electron loops, 
and, at least for s-states, for the effect of finite nuclear size. This has been done (a long 
time ago) [9l [65l [66l [67] for the hyperfine structure, but up to now not for the "muon 
Lamb shift". 

Calculations of the correction due to a vacuum polarization insert in the external 
photon line for hydrogen by Pachucki [4] and by Jentschura [57] do not agree. The 
contribution is estimated to be at least comparable to the experimental uncertainty, so 
this requires further investigation. The energy shifts for at least some of the terms in 
this contribution as calculated by Jentschura can be compared directly. Here vacuum 
polarization corrections to the Bethe logarithms will not be calculated, but vacuum 
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polarization corrections to the expectation values of W'^V and of — — can be calculated 

r dr 

quite easily, and a detailed numerical comparison for these parts is possible. 

The largest correction to the muon self-energy graphs involves the expectation value 
of V'^V. An effective charge density pvp for vacuum polarization can be derived from the 
Fourier transform of the Uehling potential, and was given by Eq. f l22l) . It is straightforward 
to calculate the expectation value of 4npvp{r). The result for hydrogen is 
(VVvp)2s = 4.4617 X 10"" /m-3 and (Wyp)2p = -8.7788 x IQ-^^ j^-3_ a check 
on this calculation, the correction to the contribution of a muon loop in the external 
photon line to the "muon Lamb shift" was calculated directly from the expectation value of 
4iTpvp{r). It gave precisely one half of the contribution normally labeled "mixed vacuum 
polarization", which would be expected. This number has been checked independently 
by other authors ([3]). The factor 2 comes from the fact that the two loops can appear 
in either order. The value of {V^Vvp)2s turns out to be equal to eypi (see Eq.[T6]) times 
the point Coulomb value. 

An additional correction due to distortion of the wave function of the 2s-state, was 
given by Ivanov et al. [BE] (see also [SS]), and corresponds to multiplying the point 
Coulomb value by evp2, which, as mentioned before, has the same numerical value. 
Strictly speaking, this correction does not correspond to a vacuum polarization insert in 
the external photon line, since it involves a second external photon, but it is numerically 
just as important as eypi- 

The main contribution to the muon Lamb shift has been given by ([3],|1]) 



m s 11 3n 

'{ZaymJ ^ 24 ^ 8^ 



Rewriting this in terms of the expectation value of V'^V, and including the contribution 
of the lowest order anomalous magnetic moment that is proportional to V'^V results in 



^Els = Tr^(VV)[[ln(— ^) + |]-lnA;o(n,/) 
37rm^ L'- ^[Zaj^rrir o 



In this work, corrections to the Bethe logarithm will not be calculated. The contributions 
that can definitely be compared with ref. [SZ] are then 



^Els,vp = Tj^(VVv^p) In ij,-^) + ^ (23) 



37rm^ 



m X 5' 



{Zaymr 6 



Calculating the energy shifts with the values calculated from Eq.[23]for (V^VVp) gave en- 
ergy shifts of 0.00191 meV for the 2s state and -0.00004 meV for the 2pi/2 state. Including 
the contribution from evp2 increases the shift for the 2s-state to 0.00486 meV for a total 
contribution of -0.00490 meV. Including the effect of finite size on the Uehling potential 
reduced the number for the 2s state by 0.00006 meV. 

In Eqs. (3.6) and (3.14) of Ref. [57] the 10/9 should be replaced by 5/6 in order to 
ensure consistency with Eq. (2.28). Doing that, and calculating the expression for the 
same pieces of Eq. (3.14) as in Eq.[2ni which is 



Za^ tZarrir^?,^ ^ r, / ni ^ 5 

2 />>') 2 
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with Vqi = 3.09 for the 2s state and Vqi =-0.023 for the 2pi/2 state gives energy shifts of 
0.004885 meV for the 2s state and -0.000036 meV for the 2pi/2 state. These numbers agree 
fairly well with those calculated here. 

The contribution to the spin-orbit shift due to the muon's anomalous magnetic mo- 
ment, with a„ = a/2TT can also be checked. The expression for — — including the Uehling 

r dr 

potential is given by Eq. f lT2|) . e2p was defined as the ratio of the value of the expectation 
value for the Uehling potential to that for the point Coulomb potential [aZ/r^), and an 
expression is given in Eq. (1131) . The numerical integration is straightforward and has been 
checked for accuracy. The result for hydrogen is = 0.000365. Thus, to obtain the energy 
shift of the 2pi/2 state, one simply has to multiply the standard value (-0.01 1713 me V) 
by £2p, giving —4.275 x 10~^meV. Evaluating Jentschura's Eq.(3.10), which neglected a 
factor rrif^/mr, gave —3.2265 x 10~^meV. Correcting for the missing mass factor gave 
-3.590 X 10-6 meV. 

For an estimate of the effect of nuclear size, an exponential charge density is used, 
mainly because all integrals can then be calculated analytically. The charge density is 
then 

o(r) = 3_e-''Ao 

Evaluation of the expectation value gives (expanded to lowest order in {aZnirro)) 
(VVc)2s = -^^^l^[l-6(aZm,ro) + 21(aZm,ro)2] 



and 

(VVc)2p = — ^ — '-{aZmrroY 

Note that aZmrVo ~ 0.00087 for muonic hydrogen, giving a value for {GaZmrTo) of 
0.0052. As was demonstrated a long time ago [67], the correction due to finite nuclear 
size is probably given correctly by 2aZmr{r)2 (here with charge density rather than with 
magnetization density). The expectation value of V'^V for the 2p-state is proportional to 
{aZmrVo)'^ 7.7 ■ 10"'', and is thus much smaller. However, one finds (to linear order in 
aZnirro) 

'1 dV\ aZ{aZmr)^ daZmrVQ 
.r~d^/ ~ 24 ^ 2 '* 

The finite size correction to the muon self-energy of the 2s-state is then approximately 
0.00699 meV for a Zemach radius of 1.086 fm for hydrogen. The total correction to the 
muon self energy contribution (which is negative) is then 0.00207 meV, which is too small 
to explain the discrepancy in the measurement of the proton radius. 

Since all of these corrections have also been calculated for deuterium and the helium 
isotopes, one can give a set of partial corrections also for these cases. Numerical values 
for the corrections eypi + £vp2 and ezem for the four cases discussed here are given in 
Table M 

Corrections to the energy shift of the 2s-state due to these effects are summarized in 
TablelTBl All energies are given in meV. The standard value is calculated with Eq. (123 p . 
and does not include all contributions, since corrections to the term involving the Bethe 
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Hydrogen 


0.00537 


-0.00762 


Deuterium 


0.00555 


-0.01918 


3He 


0.00821 


-0.03619 


^He 


0.00824 


-0.03315 



Table 15: {eypi + eyp2) and ez, 





standard value 


correction (VP) 


correction (size) 


corrected value 


Hydrogen 


-0.91689 


-0.004924 


0.00699 


-0.91480 


Deuterium 


-1.06116 


-0.00589 


0.02035 


-1.04670 


3He 


-15.5515 


-0.1277 


0.5630 


-15.1162 


^He 


-15.9455 


-0.1314 


0.5294 


-15.5475 



Table 16: total corrections to the muon self energy energy shift of the 2s-state 



logarithm are not included. The effect is significant for the helium isotopes. These 
corrections were not included in Table|2i 

As mentioned before, the spin-orbit terms in the "muon Lamb shift" should also be cor- 
rected for vacuum polarization using the correction calculated in Eq. (fT3|) for the 2p-states. 
For hydrogen, = 0.000365, for deuterium, = 0.000391, for ^He, £2^ = 0.000894, and 
for ^He, £:2p = 0.000902. The fine structure splitting would be increased by 0.00001 meV 
in hydrogen and by 0.0003 meV in ^He. 

Appendix D: Relativistic calculation of the hyperfine structure of the 2s state 

A recent paper by Carroll et al. |59] calculates the Lamb shift in muonic hydrogen rela- 
tivistically and nonperturbatively. Their calculation of the magnetic hyperfine structure 
of the 2s-state approximates the magnetization density as a delta function. However the 
magnetization density has a finite extent, and vacuum polarization has an effect on the 
magnetic vector potential. To estimate these effects relativistically, they will be calculated 
here with point Dirac wave functions. Carroll et al. do take into account perturbation of 
the wave function, which is also important. A completely correct calculation takes both 
contributions into account. 

According to Eq. (40) of ref . [T] , the hyperfine structure of an ns-state is given by 

Attko a dr r 

A^nK = , [F(F+l)-S2(s2 + l)-3/4]- / —Fn.{r)Gn.{r)- drNr%ix{rN) 

K — i/4:Zmi\f Jq r Jq 

(24) 

where and G„k are the small and large components of the wave function, n is the 
principle quantum number and n = £{£ + 1) — j(j -|- 1) — 1/4. Nuclear magnetic moments 
are given in units of nuclear magnetons, so that the proton has g = 2(1 -|- Hp). Here the 
magnetization density is normalized to 
4:71 drNr%^{rN) = 1. 

The paper by Carroll et al. [59] uses the approximation An Jq drj^rj^fi{r]\r) = 1. 
However the magnetization density has contributions from finite extent and from vacuum 
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polarization that must be taken into account. Here this will be cone perturbatively, using 
point Dirac wave functions. The wave functions are given in the book of Akhiezer and 
Berestetskii [31] . 

For an estimate of the effect of a finite extent on the magnetization density an 
exponential magnetization charge density is used, mainly because all integrals can then 
be calculated analytically. The density is then 

Here ro corresponds to the magnetic rms radius of the nucleus, and is not necessarily the 
same as that corresponding to the charge radius. 

Vacuum polarization affects the magnetic vector potential in the same manner as it 
does the electrostatic potential, and the contribution to the magnetization density is thus 
essentially the same as Eg. [2^ 



2 

TT 



q^U2iq)GM{q)joiqr) dq 



where U2{q) is defined in pj. See also Eq. The finite extent of the magnetization 
density can be taken into account when the momentum space representation is used. As 
usual, Gm is the magnetic form factor. One then finds 



Ait 



\ tq z To / 

2a r - 1)'^' M /1 o A 2m rz^ 

+ ^-^.^l + -j.(l + 2.n,r.).e-— 

(25) 

If one uses the momentum space representation, one finds for the vacuum polarization 
contribution 

47r / drNrlfPvpirN) = rcos(gr) U2{q)GM{q)dq 

Jo Jo \ q J 

The expression given in Eq.[25]can be derived directly from this form if Gm = 1- 

It is interesting to notice that if one were to define a "magnetic" potential according 
to V^Vm = 47r/i(r), then 

An drNr%n{rN) = ^'^-^ 

This permits a more transparent comparison with results obtained by other methods. It is 
easy to see that the vacuum polarization contribution to Eq.[25]is identical to the vacuum 
polarization contribution to 1'") given in Appendix A. 
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If the HFS given by Eq.[21]is calculated using An dr]\fr'j^fi{r]\f) = 1, corresponding 
to the first term in the expression for a finite magnetization density (Eq. [25|) . for the 
2s-state with j=l/2 and k = —1, the result is 



AE. 



2s 



2s 



27- 1 



(27 + 1)2 



(26) 



The factor C2S is defined as 

ag TTlr 



2s 



6mp ANiiNi + 1) 1(27 + 1) 



3/2] 



with 7 = a/1 — {aZy, Ni = + 7), and b = 2aZmr/Ni. A factor nir/m^ has 

been included to account for the fact that the magnetic moment of the muon is defined 
in terms of the free space mass, and not the reduced mass. It was pointed out in ref. pQ 
that the analysis given there does not take mass corrections to the hyperfine structure 
into account. Similar mass corrections were made by Carroll et al. [59j . 

If this is expanded in powers of aZ, and only the leading nonvanishing terms are 
retained, the result is 



c?Zg 



2s 



12m ^nip 



{2aZmr/nf [F{F + 1) - 3/2] 



Comparing this with Eq.[T3l it is easy to see that the factor multiplying [F{F + 1) — 3/2] 
is —f3/2 if Z=l and g = 2(1 + Kp). Thus, except for the missing factor 1 + a^, this is 
nearly the same as the standard perturbative result. Numerically, the value for the total 
hyperfine splitting calculated here is 22.8079 meV; this is equal to (1 + eBreit)(3 where 
P is given by Eq.[T31 as expected, since the Breit correction was defined in terms of an 
expansion of the expression for Ai?2s given in Eq.|2n]in powers of [aZ)"^, as pointed out 
in ref. [3]. Carroll et al. obtain 22.8229 meV for this contribution, which is not corrected 
for distortion of the wave functions due to vacuum polarization or finite nuclear size. The 
additional contribution due to distortion of the wave function due to vacuum polarization 
calculated here (0.0744 meV) is in fair agreement with that calculated by Carroll et al. 
(0.0747 meV). 

The magnetization density also has to be corrected for vacuum polarization. A 
calculation of Eq.|21]for the vacuum polarization correction to equation[25] gives for the 
energy shift 



AE. 



2s 



c, 



2s 



2a 
3^ 



(^2 _ 1)1/2 



-N,{N, + 2){- 



1 + 



2^2; (1 + 02)27-1 



+ 



az 



+ 



-27 -1 l + az 
L l + az^ 



) 



(27 + 1) (l + az 



2az 



az 



27 + 1 l + az 



)] 



dz 



where a = 2me/b. The resulting contribution to the total hyperfine splitting of 0.0481 meV 
agrees very well with the perturbative value of 0.0482 meV calculated nonrelativistically. 
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Including the effect of finite size on tlie vacuum polarization contribution to tlie magneti- 
zation density reduced eacli of tlie nonrelativistic perturbative contributions (modification 
of magnetization density and perturbation of the wave function) by — 0.00114 meV. 

To estimate the effect of finite extension of the magnetization density, the integral for 
the ro-dependent terms in Eq.[25]is evaluated. The resulting correction to the energy shift 
can be expanded in powers of bro. If only the leading term is retained, the result is 

AE2S = C2s{brof^-'[-N^{N^ + 2){^-^ + l + J) 

The correction to the hyperfine structure of the 2s-state due to finite extent of the magne- 
tization density is proportional to (fero)^'*'"^; it has thus been shown to be predominantly 
linear in the radius, with smaller corrections with higher powers of bro. Corrections due 
to distortion of the wave functions due to the Coulomb potential in a nonrelativistic 
perturbative calculation were historically approximated by using the Zemach radius (see 
ref. [3]). In a nonperturbative calculation, a fit to a formula that is linear plus quadratic in 
the radius would be a useful approach. In any case, the finite extent of the magnetization 
density has a dominant contribution that is linear in the radius parameter. 

The value of ro is not certain. One possibility is to take tq = rmag/V^', another is to 
use ro = 8r2em/35. The radius of the magnetization density is not necessarily the same 
as the charge radius. The magnetic radius from the Mainz experiment [17] was given as 
0.777(17) fm, which has recently been modified to 0.803(17) fm. A recent determination 
by the group at Jefferson Laboratory [TB] gives 0.867(20) fm. Numerically, one finds a 
reduction of the total hyperfine splitting of —0.111 meV if the (revised) magnetic radius 
from the Mainz experiment is used, and — 0.120 meV if the result from Jefferson 
Laboratory [16] is used. Both differ from the usual perturbative result (approximately 
-0.17meV), for reasons that are unclear. 

In any case, the correction estimated here differs significantly from the range of 
standard values obtained with a nonrelativistic, perturbative calculation, so that one 
can only say with confidence that the correction is dominated by a term linear in the 
radius parameter. 
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